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ABSTRACT 

We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way, both 
from observations and from theory under various assumption about the relation between cos- 
mic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in 
the Galactic magnetic field suggested by the Faraday rotation and polarization data, the ob- 
servations are inconsistent with local energy equipartition between cosmic rays and magnetic 
fields. Our analysis of synchrotron fluctuations suggests that the distribution of cosmic rays 
is nearly uniform at scales of order 100 pc, in contrast to that of the interstellar magnetic 
field. A conservative upper limit on the relative variations in the cosmic ray number density is 
0.2-0.4 at the scales of order 100 pc. Our results are consistent with a mild anticorrelation be- 
tween cosmic-ray and magnetic energy densities at these scales. Energy equipartition between 
cosmic rays and magnetic fields still may hold, but at scales exceeding 1 kpc. 

Key words: cosmic rays - magnetic fields - galaxies: ISM - galaxies: magnetic fields - radio 
continuum: galaxies - radio continuum: general 



1 MOTIVATION AND BACKGROUND 

The concept of energy equipartition between cosmic rays and 
ma gnetic fields and similar assumptions such as pressure eq ual- 
ity iLongaiJll994lBeck & Krausell200^ ; I Arbutina et alj|2012l) are 
often used in the analysis and interpretation of radio astronomical 
observations. This idea was originally suggested in order to esti- 
mate, from a measurement of the synchrotron brightness of a radio 
source, the magnetic field an d cosmic ray energies of the source 
as a whole (Burbidge 1956b. a). A physically attractive feature of 
the equipartition state is that it approximately minimizes the total 
energy of the radio source. 

The energy density of cosmic rays is mainly determined by 
their proton component, whereas the synchrotron intensity depends 
on the number density of relativistic electrons. Therefore, in order 
to estimate magnetic field energy, an assumption needs to be made 
about the ratio of the energy densities of the relativistic protons and 
electrons; the often ado pted value for this rati o is 100, as suggested 
by the Milky Way data dBeck & K rause 2005). This ratio is adopted 
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to be unity in applications to galaxy clusters, radio galaxies and 
active objects dCarilli & Tavlorfl2002l) . 

However, more recently this concept has been extended to 
large-scale trends in synchrotron intensity and to local energy den- 
sities at sub-kiloparsec scales in well-re s olved radio sources, such 
as spiral galaxies (e.g . iBeck et al.ll2005l:lBeckll2007l : Ichvzvll2008l : 
Tabatabaei et al. 200a lFletcher et al.l201 ll)~ Another import ant ap- 
plication of t he eq uipartition hypothesis, first suggested by Parked 
jl966|. 1 19691 , 1 19791) , is to the hydrostatic equilibrium of the inter- 
stellar gas. Here magnetic and cosmic ray pressures are assumed 
to be in a constant ratio, in practice taken to be unity. This appli- 
cation appeals to equipartition (or, more precisely, pressure equal- 
ity) at larger scales of the order of kiloparsec. The spatial relation 
between fluctuations in magnetic field and cosmic rays is crucial 
for a proposed method to measur e magnetic helicity in the I SM 
(Oppermann et al. 201 1) (see also Volegova & Stepanov 2010). 

The physical basis of the equipartition assumption remains 
elusive. Since cosmic rays are confined within a radio source by 
magnetic fields, it seems natural to expect that the two energy den- 
sities are somehow related: if the magnetic field energy density es 
is smaller than that of the cosmic rays, e cr , the cosmic rays would 
be able to 'break through' the magnetic field and escape; whereas 
a larger magnetic energy density would result in the accumula- 
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tion of cosmic rays. Thus, the system is likely to be self-regulated 
to energy equipartition, es ~ £cr- A slightly different version of 
these arguments refers to the equality of the two pressures^] giving 

However plausible one finds these arguments, it is difficult to 
substantiate them. In particular, models of cosmic ray confinement 
suggest that the cosmic ray diffusion tensor depends on the ratio 
(SB/Bo) 2 , where SB is the magnitude of magnetic field fluctua- 
tions at a scale equal to the proton gyroradiu s and Bq is the mean 
magnetic field (e.g.. iBerezinskii et alJll990l) . The magnetic field 
strength can determine the streaming velocity of cosmic rays via 
the Alfven speed, but the theory of cosmic ray propagation and 
confinement relates e cr to the intensity of cosmic ray sources rather 
than to the local magnetic field strength. Despite their uncertain 
basis, equipartition arguments remain popular as they provide 'rea- 
sonable' estimates of magnetic fields in radio sources, and also be- 
cause they often offer the only practical way to obtain such esti- 
mates. 

Equipartition between cos mic rays and magnetic fi elds can 
rarely be tested observationally. IChi & Wolfendalel d 1993b used 7- 
ray observations of the Magellanic clouds to calculate the energy 
density of cosmic rays independently of the equipartition assump- 
tion. They further calculated magnetic energy density from radio 
continuum data at a wavelength of about A12 cm. The resulting 
magnetic energy density is two orders of magni tude larger than 
that of cosmic rays, and IChi& Wolfendalel d 1993b argue that the 
discrepancy cannot be removed by assuming a proton-to-electron 
ratio for cosmic rays different from the standard value of 100 (see, 
however. Ipohll 1993b . 

An independent estimate of magnetic field strength can be ob- 
tained for synchrotron sources of high surface brightness (e.g., ac- 
tive galactic nuclei) where the relativistic plasma absorbs an ob- 
servable lower-frequency part of the radio emission (synchrotron 
self-absorption). Then magnetic field strength can be estimated 
from the frequency, the flux density and the angular size o f the syn- 
chrot r on source at the turnover f r equency (Slish 196 31 IWilliamsJ 
19631: IScheuer & Williams! Il968b . IScott & Readheadl j 1977b "and 



Read head dl994¥ concluded, from low-frequency observations of 



compact radio sources whose angular size can be determined 
from interplanetary scintillations, that there is no significant ev- 
idence of strong departures from equipartition. In the sources 
with strong synchrotron self-absorption in their sample, the to- 
tal energy is within a fac tor of 10 above the minimum energy. 
lOrienti & Dallacasa (2008) observed, using VLBI, five young, very 
compact radio sources to suggest that magnetic fields in them are 
quite close to the equipartition value. Physical conditions in spi- 
ral galaxies are quite different from those in compact, active radio 
sources, and departures from equipartition by a factor of several in 
terms of magnetic field strength would be quite significant in the 
context of spiral galaxies. 

Here we tes t the equipartition hy pothesis using another ap- 
proach (see also Step anov et all [2009). We calculate the relative 
magnitude of fluctuations in synchrotron intensity using model ran- 
dom magnetic field and cosmic ray distributions with a prescribed 
degree of cross-correlations. When the results are compared with 



observations, it becomes clear that local energy equipartition is 
implausible as it would produce stronger fluctuations of the syn- 
chrotron emissivity than are observed. Instead, the observed syn- 
chrotron fluctuations suggest weak variations in the cosmic ray 
number density or an anticorrelation between cosmic rays and mag- 
netic fields, perhaps indicative of pressure equilibrium. We con- 
clude that local energy equipartition is unlikely in spiral galaxies 
at the integral scale of the fluctuations, of order 100 pc. We discuss 
the dynamics of cosmic rays to argue in favour of equipartition at 
larger scales of order 1 kpc, comparable to the scale of the mean 
magnetic field and to the cosmic-ray diffusion scale. 



2 THE MAGNITUDE OF FLUCTUATIONS IN 
INTERSTELLAR MAGNETIC FIELDS 

Interstellar magnetic fields are turbulent, with the ratio of the ran- 
dom magnetic field to its mean component known from observa- 
tions of Faraday rotation, independently of the equipartition as- 
sumption. Denoting the standard deviation of the turbulent mag- 
netic field by of = B 2 — Bq and the mean field strength as 
Bo = I B I , where bar denotes appropriate averaging (usually vol- 
ume or line-of-sight averaging), the relative fluctuations in mag- 
netic fiel d strength in the Solar vicinity of the Milky Way is esti- 
mated as jRuzmaikin e t aljl98 S lOhno & Shibatall993l : lBeck et alj 
1996) 
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(1) 



Similar estimates result from radio observations of nearby spi- 
ral galaxies where the degree of polarization of the integrated emis- 
sion at 4.8 GHz is a few per cent, with a range p ~ 0.01-0.18 
dStil et alj r2009). These data are affected by beam depolarization, 
so they only give upper limits for 8 2 . More typical values of the 
fractional polarisation in spiral galaxies are p — 0.01-0.05 within 
spiral arms and 0. 1 on average. The degree of polarization at short 
wa velengths where Faraday rotation is negligible can be estimated 
as teurnll 1 9661: ISokoloff et al Jl 1 998b 
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where Box is the strength of the large-scale magnetic field in the 
sky plane, po w 0.75, and the random magnetic field is assumed to 
be isotropic, b\ = '4 a 2 . This yields 
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in a good agreement with the Milky Way data obtained from Fara- 
day rotation measures. Forp = 0.05-0.1, we obtain 5 2 ~ 10-20. 

It is important to note that Eq. ((3) has been obtained assuming 
that the cosmic ray number density n CI is uniform, so that all the 
beam depolarization is attributed solely to the fluct uations in mag- 
netic field. Under local equipartition, n CI oc B 2 , ISokoloff et all 
( 1998, their Eq. (28)) calculated the degree of polarization at short 
wavelengths to be 



It is useful to carefully distinguish between what can be called 'pressure 
equality' and 'pressure equilibrium': the former refers to the case where 
magnetic fields and cosmic rays have equal pressures locally, whereas the 
latter describes the situation where the sum of the two (or more) pressure 
contributions does not vary in space. 
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As might be expected, this expression leads to a smaller 8b for a 
given p/po than Eq. $2%: 
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so that 
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trie radius), I and v are the correlation length[f] and r.m.s. speed of 
the interstellar turbulence, so that l/v is the lifetime of a turbulent 
eddy. This leads to 



Since local equipartition between cosmic rays and magnetic fields 
maximizes beam depolarization, this is clearly a lower estimate of 

5 b . 

The above estimates apply to statistically isotropic random 
magnetic fields. However, the random part of the interstellar mag- 
netic field can be expected to be anisotropic at scales of or- 
der 100 pc, e.g., due to shearing by the galactic differential ro- 
tation, streaming motions and large-scale compression. Note that 
the anisotropy of MHD t urbulence resulting from the nature of the 
spectral energy cascade jLithwick et alj|2007t IcTaltier et alj|2000l 
and references therein) is important only at much smaller scales. 
Synchrotron emi ssion arising in an anisotropic ran dom magnetic 
field is polarized (Laing 1980; Sokolof f et all 19981) and the result- 
ing net polarization, from the combined random and mean field, 
can be either stronger or weaker than in the case of an isotropic 
random field depending on the sense of anisotropy relative to the 
orientation of the mean magnetic field. 

The case of M33 provides a suitable illustration of the refine- 
ments r equired if the anisotropy of the random magnetic field is sig- 
nificant. |Tabatateej_rt_alJ ([2008J) obtained fractional polarization of 
about 0.1 at A3. 6 cm. Using Eq. {3), this yields <5 2 ~ 10, whereas 
Eq. (01 leads to 5' 2 ~ 4, consistent with their equipartition estimates 
at — 6 /iG and Bo — 2.5 /iG. However, their analysis of Faraday 
rotation between A3. 6, 6.2 and 20 cm suggests a weaker regular 
magnetic field, Bo — 1 /iG, leading to 5 2 ~ 40 if ab — 6 /iG. 

The latter estimate for Bo is more reliable since using the de- 
gree of polarization can lead to an un derestimated 5b becau se of 
anisotropy of the random magnetic field. Sokoloff et al.l d 19981 their 
Eq. (19)) have shown that the degree of polarization at short wave- 
lengths in a partially ordered, anisotropic magnetic field is given 

by 



■ = Po 
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where we have assumed, for the sake of simplicity, that B y — Bq 
and B x = 0, defined 5 2 = <r 2 /B 2 and likewise for 5b x , and intro- 
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/a by (< 1) as a measure of the anisotropy of bj 



This approximation is relevant to spiral galaxies where the mean 
magnetic field is predominantly azimuthal (nearly aligned with the 
y-axis of the local reference frame used here) and the anisotropy 
in the random magnetic field is produced by the rotational shear, 
<Jb y > Cbx- For Sby 3> 1, this yields, for p — 0.1, 
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Thus, a rather weak anisotropy of the random magnetic field can 
produce p ~ 0.1 and this allows us to reconcile the different es- 
timates of S b obtained from the degree of polarization and Fara- 
day rotation in M33. The required anisotropy can be readily pro- 
duced by the galactic differential rotation. Shearing of an initially 
isotropic random magnetic field leads, within one eddy turnover 
time, to an increase of its azimuthal component of 
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where fi(r) is the angular velocity of the galactic rotation (with the 
rotational velocity along the local y-direction and r the galactocen- 
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where the last equality is based on the estimate rdfl/dr ~ 
— Vo/-Ro> with Vo = 107kms _1 and Ro = 8kpc being the pa- 
rameters of Brandt's approximation to the rotation curve of M33 
dRogstadetaIlll976l) . With I = 0.1 kpc and v = lOkms" 1 , 
we obtain a 2 , ~ 0.8, in a perfect agreement wi th the degree of 
anisot ropy required to explain the observations of iTabatabaei et all 
( l2008h . 

To conclude, a typical value of the relative strength of the ran- 
dom magnetic field in spiral galaxies is, at least, 
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(5) 



This estimate refers to the correlation scale of interstellar turbu- 
lence, I ~ 50-100 pc. Higher values, <5 2 ~ 40 are perhaps more 
plausible, but our results are not very sensitive to this difference 
(see Fig.[T]and Section[3jl. 



3 SYNCHROTRON INTENSITY FLUCTUATIONS 

The ratio of the fiuctuating-to-mean synchrotron intensity in the in- 
terstellar medium (ISM) is sensitive to the relative distributions of 
cosmic ray electrons and magnetic fields and hence to the extent 
that energy equipartition may hold locally: the synchrotron emis- 
sion will fluctuate strongly if equipartition holds pointwise, i.e., if 
the number density of cosmic ray electrons is increased where the 
local magnetic field is stronger. (We assume that cosmic ray elec- 
trons and heavier particles are similarly distributed - see Section[8] 
for the justification.) We estimate the area-averaged synchrotron in- 
tensity Iq = (J) g in the plane of the sky and the standard deviation 



of the intensity fluctuations, cri = y (J 2 ) s — (I) s , from the sta- 
tistically homogeneous, diffuse ISM. Assuming for simplicity that 
the synchrotron spectral index is equal to — 1, the total intensity of 
synchrotron emission is given by 



e ds oc / n cr B x ds, 
l J L 



(6) 



where e is the synchrotron emissivity, n cr is the number density of 
cosmic ray electrons, B± is magnetic field in the plane of the sky 
and integration is carried along the line of sight s over the path- 
length L. 

The relative fluctuations in the synchrotron intensity in the 
case of detailed (local, or pointwise) equipartition between cosmic 
rays and magnetic fields, n cr oc B 2 can be estimated as (see Ap- 
pendix |Aj 



oj_ _ <5 h (54 + 295<S 2 + 404<^ + 101£f ) 1/2 
Jo ~ N 1 / 2 (3 + + 5<5 b 4 ) 



(7) 



where is the number of correlation cells of the synchrotron fluc- 
tuations within the telescope beam, and each vector component of 



2 The correlation length is also known as the integral scale and is defined as 
the integral of the variance-normalized autocorrelation function of a random 
variable. The typical linear size, or diameter, of a turbulent cell is 21. 
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Figure 1. The relative fluctuations in synchrotron intensity 5j for spe- 
cial cases: complete correlation of cosmic rays and magnetic fields (solid), 
n CT = const with fluctuating b (dashed) and uncorrected fluctuations in 
cosmic rays and magnetic field, with 5 n = 1 (dotted) and 8 n = 0.5 (dot- 
dashed). The corresponding analytical expressions can be obtained from the 
formulae presented in Appendix [Aland Section fJ] We note that the curves 
rapidly approach the horizontal asymptote, and the approximation <5j, — !• oo 
is reasonably accurate for <5{, > 1.5-2. 



the magnetic field is assumed to be an independent Gaussian ran- 
dom variable of the same variance (i.e. an isotropic random mag- 
netic field), and we recall that 5b = <Jb/Bo. As St increases, Si 
rapidly approaches the asymptotic value 

Si « 2N~ 1/2 for 5 2 b > 1. 

It is useful to note similar expressions for Si obtained under 
different assumptions about the correlation between cosmic rays 
and magnetic field. If cosmic ray fluctuations are statistically inde- 
pendent of those in magnetic field, with the mean value no = n cr 
and standard deviation <r n , we show in Appendix lAl that 



[5 2 n + 5 2 (2 + 5 2 )(l + 25 2 )} 1 / 2 
N 1 / 2 (1 + 51) 



(8) 



where 



On = 

n 



is the magnitude of the relative fluctuations in cosmic ray number 
density. In particular, for 8 n ->0we obtain an asymptotic form for 
a homogeneous distribution of cosmic rays: 



Si = 



5 b (2 + 8tf' 2 



(9) 



iVV2 (i + p b ) ■ 

Thus, 8i ~ TV~ 1/2 for 5 2 > 1 and n cr = const. 

For the correlation length of the synchrotron fluctuations l e — 
50 pc and the pathlength L — lkpc, we obtain TV = L/2l e ~ 
10 for a beam narrower than the size of the correlation cell. With 
Eq. J5j, the relative magnitude of synchrotron fluctuations expected 
under detailed equipartition follows from Eq. (|7} as 

N x ^ 2 8i w 2.0 (local equipartition). 



(10) 



Equation yields, for 5„ = 0.5, 

7V 1//2 5_f « 1.2 (uncorrected fluctuations), 
and Eq. $9^ leads to 

N 1/2 Si « 1.0 (n c 



; const) 



As might be expected, detailed equipartition between cosmic rays 
and magnetic fields leads to the strongest synchrotron fluctuations 
for a given 5b and TV. Figure [T] shows the dependencies defined by 



Eqs 0, ^ and ((9}. We note that the dependence of Si on 5b is 
quite weak as long as 5 2 > 3. 



4 SYNCHROTRON FLUCTUATIONS IN THE MILKY 
WAY AND M33 

In this section we estimate the relative level of synchrotron fluctu- 
ations from observations of the Milky Way and the spiral galaxy 
M33. An ideal data set for this analysis should (i) resolve the fluc- 
tuations at their largest scale, (ii) only include emission from the 
ISM and not from discrete sources such as AGN and stars, (iii) not 
be dominated by structures that are large and bright due to their 
proximity, such as the North Polar Spur, (iv) be free of systematic 
trends such as arm-interarm variations or vertical stratification. The 
data should allow the ratio 5i = <ti/ Jo to be calculated separately 
in arm and inter-arm regions or at low and high latitudes as Jo dif- 
fers between these regions. Regarding item (i) above, we note that 
a turbulent cell of 100 pc in size subtends the angle of about 6° at 
a distance of 1 kpc. Furthermore, most useful for our purposes are 
long wavelengths where the contribution of thermal radio emission 
is minimum. Unfortunately, ideal data satisfying all these criteria 
do not exist; we therefore use several radio maps, where each map 
possesses a few of the desirable properties listed above and collec- 
tively they have them all. 



4.1 The data 

4.1.1 The 408 MHz all-sky survey 

The survey of lHaslam etai] dl982h covers the entire sky at a res- 
olution of 51' and with an estimated noise level of about 0.67 K. 
Synchrotron radiation is the dominant contribution to emission at 
the survey's wavelength of A74 cm. The brightest structures in the 
map shown in Fig.|2ji are the Galactic plane and several arcs due to 
nearby objects, especially the North Polar Spur. 

4.1.2 The 408 MHz all-sky survey, without discrete sources 

lLa Porta et al. 1 20081) rem oved the strongest discrete sources from 
the data of Hasl am et "all Jl982h using a two-dimensional Gaussian 
filter. We compared the results obtained from the original 408 MHz 
survey with those from this map to show that the effect of point 
sources on our results is negligible. 



4.1.3 The 22 MHz part-sky survey 

iRoger et alj dl999l) produced a map, shown in Fig. [3^, of about 
73% of the sky at A13.6m, between declinations —28° and +80° 
at a resolution of approximately 1° x 2° and an estimated noise 
level of 5 kK. The emission is all synchrotron radiation, but H II 
regions in the Galactic plane absorb some background emission 
at this low frequency. However, we are most interested in regions 
away from the Galactic plane, so our conclusions are not affected 
by the absorptio n in the H II region s. The brightest point sources 
were removed bv lRoger et al.l d 19991) as they produced strong side- 
lobe contamination in the maps: this accounts for the four empty 
rectangles in Fig.|3ji. 

We expect that results useful for our purposes arise at the an- 
gular scale of about 6° in all three Milky Way maps (i.e., the angu- 
lar size of a turbulent cell at a 1 kpc distance), whereas larger scales 
isolate regular spatial variations of the radio intensity. 
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4. 1.4 The 1.4 GHz map ofM33 

The nearby, moderately inclined, spira l galaxy M33 was observed 
at A21 cm bv lTabatabaei et al] ( 1200 7a'). using the VLA and Effels- 
berg 100 m telescopes, at a resolution of 51", or about 200 pc at 
the distance to M33 of 840 kpc. The noise level is estimated to 
be 0.07 mjy/beam. The resolution is sufficient to resolve arm and 
inter-arm regions, but is at the top end of the expected scale of ran- 
dom fluctuations due to turbulence. The beam area includes a few 
(nominally, four) correlation cells of the synchrotron fluctuations. 
The emission is a mixture of thermal and synchrotron radiation. 
The overall thermal fraction is estimated to be 18 % but it is strongly 
enhanc ed in large Hll regions and spiral arms (Tabatab aei et al.l 
2007b) whereas the synchrotron emission comes from the whole 
disc. The radio map used here is shown in Fig.|4^. The spiral pattern 
in notably weak in total synchrotron intensity, so the map appears 
almost featureless. This makes this galaxy especially well suited for 
our analysis since we are interested in quasi-homogeneous random 
fluctuations of the radio intensity. Nevertheless, systematic trends 
are noticeable in this map and we discuss their removal in Sec- 
tion l4273l 

4.2 Statistical parameters of the synchrotron fluctuations 

For the three Milky Way data sets of Sections [4. 1 . 1444. 1 .31 we cal- 
culated the mean Io and standard deviation 07 of the synchrotron 
intensity / at each point in the map. In each case, the data were 
smoothed with a Gaussian of an angular width a, resulting in the 
local mean intensity at the scale a, which we denote Io a : 

Ioa = S~ 1 J J I(l',b')exp(-8 2 /a 2 )cosbdl'db', (11) 

where integration extends over the data area, r = (1, 1, b) is the 
position vector on the unit sky sphere, with I and b the Galactic lon- 
gitude and latitude (confusion with the small-scale magnetic field, 
denoted here b, should be avoided), 8 = arccos(r • r') is the an- 
gular separation between r and r', 

S a (l,b) = jj exp(~0 2 /a 2 )cosbdl'db', 

is the averaging area, and the integration extends over the whole 
area of the sky available in a given survey. The standard deviation 
a ia of radio intensity at a given position (I, b) at a given scale a is 
calculated as 

a 2 Ia (l,b) = (l 2 ) a -{I)l, 

where angular brackets denote averaging defined as in Eq. (lilt . 

In the case of M33, we selected nine areas which avoid the 
brightest H II regions and whose radio continuum emission is thus 
likely to be dominated by synchrotron radiation. Each area encom- 
passes several beams and 81 was calculated in each, using the mean 
value and the standard deviation of / among all the pixels in the 
field obtained after removing regular trends. 

4.2.1 The 408 MHz survey 

To reduce the influence of the Galactic disc, where the compli- 
cated structure in the radio maps is mainly due to systematic arm- 
interarm variations and localized radio sources such as supernova 
remnants, the original intensity data were truncated at 52 K (1% of 
the maximum and 167% of the r.m.s. value of /). The resulting sky 
distribution of the relative radio intensity fluctuations oi a /I<3a is 
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Figure 2. (a): The 4 08 MHz all-sky map of the total synchrotron intensity 
IHaslam et alj|l982fl . with the Galactic disc area (I > 52 K) blanked out. 
The lower panels show the magnitude of the relative fluctuations of the 
synchrotron intensity, Sj = aj/Io, at various scales: (b) a = 30°, (c) 
a = 15° and (d) a = 7°. The latter scale is about the angular size of a 
turbulent cell (2l E = 100 pc) seen at a distance 1 kpc. 



shown in Fig. |2]for a selection of averaging scales, a — 30°, 15° 
and 7° (a corresponds to the radius rather than the diameter of the 
region). 

While panels (b) and (c) in Fig. [2] reflect mainly systematic 
trends in radio intensity, we expect that panel (d) is dominated by 
the turbulent fluctuations. In particular, the a = 7° map shows a 
much weaker variation with Galactic latitude than those at larger 
scales, for |&| > 20°. We note that the correlation scale of the syn- 
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Figure 3. As in Fig.|5]but for the 22 MHz survey jRoger et all 19991) . 
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Figure 4. The 1.4 GHz radio map of M33 iTabatabaei et al.ll2007dl with 
the rectangular fields used in our analysis shown. For orientation, we note 
that the size of Field 3 is 7' X 7' equivalent to 1.6 X 1.9 kpc 2 . 



Table 1. Relative radio intensity fluctuations oj/io in M33 with systematic 
trends of various orders subtracted. The first column gives the field number 
as specified in Fig. 0] and the next three columns show the relative fluc- 
tuations of radio intensity, with Cj" 1 ' the standard deviation of / within 
a given field, obtained with a trend of order m and zero mean value sub- 
tracted: m = corresponds to the original data, m = 1, to a linear trend, 
and m = 2, to a quadratic trend in the angular coordinates. The last column 
shows the mean value of the radio intensity in each field. 



chrotr on fluctuations obtained by IPagkesaman skii & Shutenkovl 
dl987t) is 8°. 

Contours outside the Galactic disc in Fig.|2jl, |6| > 20° give 
= 0.1-0.2. (12) 

lo 

Results obtained from the 408 MHz map with point sources 
subtracted differ insignificantly from those obtained using the orig- 
inal map. 



Field No. 


4 0) //o 




4 2) //o 


Iq [/i.Jy/beam] 


1 


0.30 


0.18 


0.14 


677 


2 


0.34 


0.22 


0.19 


545 


3 


0.31 


0.15 


0.11 


609 


4 


0.28 


0.18 


0.17 


606 


5 


0.38 


0.17 


0.13 


672 


6 


0.35 


0.15 


0.12 


795 


7 


0.30 


0.16 


0.14 


565 


8 


0.32 


0.07 


0.07 


902 


9 


0.38 


0.14 


0.11 


810 


Mean 


0.33 


0.16 


0.13 


687 



4.2.2 The 22 MHz partial sky survey 

Contours of the relative intensity fluctuations obtained from the 
22 MHz map are shown in Fig. [3] As in Fig. [2] averaging over 
scales a — 30° and 15° reveals the large-scale structure clearly 
visible in the original data. However, results at a = 7° show much 
less of such structure, and the statistically homogeneous part of 
the sky in this panel includes the same contours of 0. 1 and 0.2 as 
in Fig. [2] with the value of 0.3 confined to the bright ridges seen 
in Fig. [3^. Thus, the 22 MHz data are in a good agreement with 
the values for Si obtained from the 408 MHz data. This suggests a 
weak frequency dependence of the relative synchrotron fluctuations 



4.2.3 M33 

Our analysis for the Milky Way has a potential difficulty that long 
lines of sight might make it impossible to separate the contributions 
to oi from small-scale (random) and large-scale (systematic) varia- 
tions of the synchrotron emissivity. Therefore, we consider also the 
nearby galaxy M33 seen nearly face-on (inclination 56°). To avoid 
excessive contribution from large-scale variations due to the spi- 
ral pattern and the radial gradient in I, we selected relatively small 
fields shown in Fig. [4] selecting areas free of strong star formation 
in the outer regions of the galaxy disc. The areas of the rectangular 
fields chosen range from 1.9 x 1.2 kpc 2 to 1.9 x 4.7 kpc 2 . 
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4.3 Comparison with earlier results and summary 




Synchrotron observations in the Milky Way focus primarily on 
the spectrum of the fluctuations while their magnitu de has at- 
tracted surprising ly little attention. iMills & Sleel jl957b observed 
fluctuations of the Galactic radio background near the Galactic 
south pole at A = 3.5 m, with the resolution of 50°, to obtain 



Figure 5. The relative fluctuations in synchrotron intensity Sj as a func- 
tion of Nw > the number of synchrotron correlation cells across the beam 
area, with Si normalized by So, the relative fluctuations obtained along 

a single line of sight, i.e., for Nw = 0. The asymptotic dependence 

— 1/2 

Si = 2SqN w 1 emerges only for Nw > 20-30. The calculation as- 
sumed a Gaussian beam, with the full width at half-maximum (FWHM) 
taken for the beam diameter. 



The fields are big enough to make gradients in the mean quan- 
tities significant; in particular, the non-t hermal disc of M33 has a 
strong radial gradient in radio intensity (Tabata baei et alj[2007bh . 
so we subtract regular trends from the values of /. We fitted first- 
or second-order polynomials to 1(1, b) in I and b in each field and 
calculated 8i after subtracting the trends with zero mean value from 
the original data. Results are shown in Table [TJ with cr^™' denot- 
ing the standard deviation of the radio intensity obtained upon the 
subtraction of a polynomial of order m in the angular coordinates. 
The mean value of synchrotron intensity Io was calculated for each 
field. 

We note that crj ^ (the standard deviation of I in the original 
data) is noticeably larger than <r} and 07 . Thus, the large-scale 
trends contribute significantly to the intensity variations. On the 
other hand, and <rj 2 ' have rather similar magnitudes of order 
0.15 (and > cr| 2 ' as expected), so that they can be adopted as 
an estimate of 07 corrected for the large-scale trends. We use the 
value of /Io averaged over the nine fields explored as the best 
estimate for Si . 

However, W w 100 pc (half-width at half maximum of the 
Gaussian beam) in the observations of M33 used here, and the 
beam area encompasses Nw ~ 4 correlation cells. For compari- 
son with the Milky Way data, this estimate has to be reduced to the 
virtually perfect resolution available in the Milky Way. For such a 
small value of Nw as in the high-resolution observations of M33, 
the dependence of Si differs significantly from its asymptotic form 
Si oc N w x l 2 . Figure [5] shows the dependence of Si on Nw ob- 
tained for a model synchrotron-emitting system described in detail 
in Section [5] Si only weakly depends on Nw for small Nw- Re- 
duced to a single line of sight, the synchrotron fluctuations in M33 
then correspond to Si — 0.2 if the synchrotron correlation length 
is l e — 50 pc (i.e., Nw = 4). We also recall that the typical path- 
length through M33 is twice that through the Milky Way, and we 
adopt L — 2 kpc for this galaxy. Then the value of Si in M33, re- 
duced to L — 1 kpc for compatibility with the Milky Way data, is 
further V2 times larger: 

Si ~ 0.3. 



to S i = 0.12 (see alsolGetmantsevlll959h. 

Dagkesamanskii & Shutenko vl dl987h used observations at 
102.5 MHz (A2.92m) near the North Galactic Pole, where the 
Galactic radio emission is minimum, to determine the synchrotron 
autocorrelation function and its anisotropy arising from the large- 
scale magnetic field. They obtain Si « 0.07 but note that this 
estimate should be doubled if the isotropic extragalactic back- 
gound (half the total flux) is to be subtracted, to yield Si ~ 0.14. 
Band av et alj dl991ij) argue that only 17% of the total flux is of 
extragalactic origin, and then Si ~ 0.08 at 102.5 MHz. 

The autocorrelation function of the brightness temperature 
flu ctuations of th e Galact ic radio background was determined also 
by iBandav et al. I <1991al l b) who used observations at 408 and 
1420 MHz, smoothed to a resolution of about FWHM = 5° . They 
observed a 'quiet' region with reduced fl uctuations, 30° < DEC < 
50°, 180° < RA < 250°, identified bv lBridld ( fl967h as an inter- 
arm region, since they were interested in the cosmic microwave 
background fluctuations. For the Galactic synchrotron radiation, 
which dominates at these frequencies, they obtain 81 ~ 0.05 at 
408 MHz and 0.08 at 1420 MHz. 

These estimates are somewhat lo wer than those obtained in 
Section [4] The value of Si obtained by IBandav et al.l can be lower 
due to their selection of a region with weaker synchrotron fluctua- 
tions. 

The relative fluctuations in radio intensity are remarkably sim- 
ilar in all the Milky Way maps and in all fields in M33 considered, 
with 

Si = ai/Io « 0.1-0.3, 

when reduced to the common pathlength L = 1 kpc (see Sec- 
tion |6j2] for further discussion). The lower values are more plau- 
sible. We believe that these estimates are not significantly affected 
by either large-scale trends in the radio intensity or by discrete radio 
sources or by thermal emission. Even if these effects still contribute 
to our estimate, it provides an upper estimate of the fluctuations 
in synchrotron intensity arising in the interstellar medium of the 
Milky Way and M33. 



5 MAGNETIC FIELD AND COSMIC RAY MODEL 

In order to interpret the results of our analysis of observations in 
Section [4] we use a model of a partially ordered random magnetic 
fields and an associated distribution of cosmic rays, assuming var- 
ious degrees of correlation (or anticorrelation) between these two 
components. We then calculate the synchrotron intensity from a 
computational box of appropriate size and compare the result with 
the observational constraints described above in order to establish 
the degree of correlation or anticorrelation compatible with obser- 
vations. 

The model magnetic field is described as a sum of a mean and 
random component, Bo and b, respectively: 

B = Bo + b, 

with volume averages (6) = 0, (B) = Bo and (B 2 ) = Bq ± + 
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+ a 2 , where B 



and B ij are the mean field components in 



the plane of the sky and parallel to the line of sight, respectively. 



2.0 - 



5.1 A model of a partially ordered magnetic field 

To prescribe a quasi-random magnetic field b with vanishing mean 
value in a periodic box, we use a Fourier expansion in modes with 
randomly chosen directions of wave vectors fe, but with amplitudes 
adjusted to reproduce any desired energy spectrum: 



b(x) = 



1 



(2tt) 3 / 2 



b(k)e ik 



(13) 



where b is the Fourier transform of b; the physical field is repre- 
sented by the real part of this complex vector. The corresponding 
magnetic energy spectrum is given by 



M(k) = 



|S(fe')| 2 d 3 fc', 



(14) 



where the integral is taken over the spherical surface of a radius 
k in the fc-space. In the isotropic case, M(k) = 4ivk 2 \b(k)\ 2 . In 
order to ensure periodicity within a computational box of size L, as 
required for the discrete Fourier transformation, the components of 
the wave vectors are restricted to be integer multiples of ^[L. 

A solenoidal vector field b, i.e., that having V • b = k ■ b(k) = 
0, is then obtained from k and M(k) as 



b(k) = 



k x X 
\k x X\ 



k-^M{k), 



where X is a complex vector chosen at random, to ensure that the 
Fourier modes have random phases. We consider a magnetic energy 
spectrum represented by two power-law ranges, 



M(k) = M 



(fc/fco)*° 
(k/k )- s 



for k < ko, 
for k ko, 



with so > and si > 0, where ko is the energy-range wavenum- 
ber. We use Si = 5/3 as in Kolmogorov's spectrum and so = 2 
(see lChristensson et alj|200lh . The standard deviation of the mag- 
netic field is given by 

so + Si 



M(k)dk = M k , 
o (so + l)(si 



1) 



for si > 1. The correlation length h of the resulting magnetic 
field (analogous to the radius of a correlation cell) differs from its 
dominant half- wavelength ^Ao = 7r/fco for any finite values of so 
and si: 



7T dk _ ji_ 

2 J °° M(k) dk ~ 2k 
3tt 



1 - 



si 



10 V 



(15) 



where the last equality follows for s = 2 and Si =5/3. 

The resulting solenoidal vector field is then added to a uni- 
form component Bo to produce a partially ordered magnetic 
field with controlled fluctuation level Si, and energy spectrum 
M(k). This approach has been used to generate synth e tic po - 
larization maps of th e turbu l ent ISM by IStepanov et all 1 200 8l): 
IVolegova & Stepanovl d2010h : lArshakian etalJfeOllh: Moss et al. 
2012 1 ) . Si milar constru c tions were used by iGiacalone & Jokipiil 
19991) and ICasse et alj fe002l) in their mode lling of cosmic ray 



propag ation in random magnetic fields, by Malik & Vas silicosl 
( 1999) for modelling turbulent flows, and by IWilkin eFaT] d2007l) 
to study dynamo action in chaotic flows. 




0.0 



-0.5 



0.0 

c 



Figure 6. The relative fluctuations in synchrotron intensity Sj as obtained 
from analytical formulae Eqs. flA3l and l lA4t ) (thicker curves) and from nu- 
merical calculations with the condition n cr > enforced (the correspond- 
ing thinner curves), for: <5 n = 0.2 (solid), 8 n = 0.4 (dashed) and 8 n = 0.8 
(dash-dotted). The analytic formulae become inapplicable for C < and 
S n — 1 (see the text). The magnetic field is purely random, Bo = and 
N = 10. 



5.2 Cosmic ray distribution partially correlated with 
magnetic field 

Now we introduce a random distribution of cosmic rays which has a 
prescribed cross-correlation coefficient with the total magnetic en- 
ergy density B 2 /8iv for the magnetic field described in Section l5~71 
We first generate another, independent realization of the magnetic 
field, B' using the same procedure as above, but with different ran- 
domly selected vectors X and B' = 0. By construction, B 2 and 
B' 2 are uncorrected scalar fields, 



where 



C{u,v) = 



C{B 2 ,B' 2 ) = 0, 

uv) — (u) (v) 



(16) 



u u a v 

is the cross-correlation coefficient of the random variables u and 
v, with C — 1 if u and v are perfectly correlated and C = — 1 
for perfectly anticorrelated u and v. The number density of cosmic 
rays is then obtained as 

n cr = 8 n [a (B 2 - (B 2 )) + (b' 2 - (s' 2 ))] + 1, (17) 

where a = C/a B 2 and /? = y/l — C 2 /a B ,2 (see Appendix lAlfor 
details). Here 



"2„2 

_ 3 



oi(2B$ + oi) 



2^1/2 



is the standard deviation of B 2 , and a B ,2 is that of B' 2 . Equa- 
tion J17t ensures that 



(n) 



n 



= S„ 



in other words, the number density of cosmic rays n cr is measured 
in the units of no. 

This model is flexible enough to admit control over the spec- 
trum of cosmic-ray density fluctuations by a suitable choice of the 
spectrum of B' . However, we do not use this opportunity, and B 
and B' always have identical spectra. 

We will now verify, by direct calculation, that C(B 2 ,n CT ) w 
C. The reason for the approximate equality is first explained. 

A shortcoming of this cosmic-ray model is that n cr can be 
negative at some positions (especially when C < and hence 
a < 0) because, at some positions and in some realizations, B 2 can 
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be arbitrarily large (as a Gaussian random variable squared). This 
deficiency could be corrected by selecting a more realistic proba- 
bility distribution for b (e.g., a truncated Gaussian) but we do not 
feel that this would lead to any additional insight. In the numeri- 
cal calculations described below, we truncate the negative values 
of n cr replacing them by zero. This, however, makes it impossible 
to achieve exact anticorrelation between cosmic rays and magnetic 
fields, so that C (B 2 , n cr ) > — 1. As shown in Fig. [6] this constraint 
is not restrictive for S n < 1, where C — — 1 is achievable with high 
accuracy; however, for 5 n — 0.8 the typical smallest value of C is 
about —0.8. Enforcing n cr = at positions where B 2 is excep- 
tionally large can lead to a reduction in the values of ai obtained. 

In analytical calculations, we restrict ourselves to the cases 
with 5 n < 1 to reduce the extent of the problem (even if not to 
resolve it completely). For example, Eqs. jA3\ and dA4t yield for 
Bo — > (i.e., S b — > oo) 



Si = 



[9 + 6S n (llC 2 S n + 3V6C + 33 n ) 



1/2 



(18) 



This dependence of Si on the cross-correlation coefficient C is 
shown with thicker curves in Fig. [fj] for various values of the rel- 
ative fluctuations in cosmic rays, 5 n . Thinner curves show similar 
results obtained from a numerical calculation where n cr > is 
enforced. It is clear from Fig. |6]that these analytical results are ac- 
curate for C > 0. However, for C < 0, they are useful only if the 
fluctuations in the cosmic ray number density are relatively weak: 
C > -0.1 for S n < 0.8, C > -0.5 for S„ < 0.4, and almost 
any value of C for S n < 0.2. We only use these analytical results 
for illustrative purposes, whereas all our conclusions are based on 
numerical results where n cr > at all positions. 



6 RESULTS 

6.1 Synthetic radio maps 

Each component of the magnetic field described by Eq. dl3t is the 
sum of a large number of independent contributions from different 
wave numbers. By virtue of the central limit theorem, each compo- 
nent of the resulting magnetic filed is well approximated by a Gaus- 
sian random variable. Then the mean synchrotron intensity and its 
standard deviation over N correlation cells can be expressed, using 
Eq. {6]l, in terms of Bo, 0b, S n and C. Explicit analytic expressions 
for 7o and 07 can be found in Appendix [A] and we illustrate these 
results in Fig. [6] As might be expected, the relative level of the syn- 
chrotron intensity fluctuations increases with the cross-correlation 
coefficient between B 2 and n cr . 

Since analytical results are of limited relevance for C < 0, 
we performed numerical calculations of the synchrotron inten- 
sity where the cosmic-ray number density is truncated to be non- 
negative (i.e. n cr = wherever the model defined by Eq. ( 117) 
returns a negative value). The model has four free parameters: 

(i) the relative level of magnetic field fluctuations St = at /Bo, 

(ii) the relative level of cosmic-ray number density fluctuations 

Sn — <Jn/no, 

(iii) the cross-correlation coefficient between magnetic field and 
cosmic rays C, and 

(iv) the dominant energy wave number of the turbulent magnetic 
field fco- 

We do not vary the spectral index of magnetic field and cosmic rays 
as these parameters are of secondary importance in this context. 




Figure 7. The isocontours of the numerical factor G in Eq. 1201 for 05 = 1 
(dashed lines and labels in squares) and u\y — > 00 (solid lines and labels in 
ovals) shown in the (<5 n , C)-plane. 



The value of fco controls the correlation lengths of magnetic 
field (Eq. J15b). cosmic rays and synchrotron emissivity, and hence 
the number of the correlation cells of synchrotron fluctuations in 
the telescope beam N, which in turn affects the magnitude of syn- 
chrotron fluctuation as Si oc N- 1/2 . Since N can vary widely be- 
tween different lines of sight in the Milky Way and between galax- 
ies with different inclination angles, we present our results in terms 
of A" 1 / 2 ^ for both the observations and the model. 



6.2 A relation between the correlation lengths of the 
synchrotron emissivity and magnetic field? 

In the case of infinitely narrow beam, the number of synchrotron 
correlation cells traversed by the emission is just the ratio N = 
L/(2l e ), where L is the pathlength through the synchrotron source 
and l e is the correlation length of the fluctuations in the synchrotron 
emissivity. For a finite beam width W, this is the number of correla- 
tion cells within the beam cylinder, N ~ (3/ 16) LW 2 / , assum- 
ing a circular beam and spherical correlation cells. Unlike the corre- 
lation lengths of physical parameters such as the magnetic field, ve- 
locity or density fluctuations, the correlation length of the intensity 
(or emissivity) variations cannot be deduced independently (e.g., 
from the nature of the turbulence driving), but has to be calculated 
from the statistical parameters of the physical variables or from ob- 
servations. 

To illustrate the difficulties arising, consider the autocorrela- 
tion function of b\ as an example. If V(x) is a stationary Gaussian 
random function, with vanishing mean value and the autocorrela- 
tion function K v (r) = V(x)V(x + r), the autocorrelation func- 
tion of W (x) = V 2 (x) is given by K w (r) = 2[K v (r)} 2 (s 



e.g. i)13 in Sveshnikov 1966). Assuming that each component bi 
of the random magnetic field is a Gaussian random variable, with 
the autocorrelation function , we have K h .i = IK 2 . . Assum- 



ing statistical isotropy of b, Kt x (r) 
cross-correlations between b x and 6, 



= Kb y (r), and neglecting any 
, we obtain the autocorrelation 
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function of b\ : 

K b 2 ± (r)=4Kl(r). 

Then relation between the correlation scales of bi and b\ , denoted 
l bi and l b 2 , respectively, depends on the form of the autocorrela- 
tion function of magnetic field: for Kb i = \o1 exp(— r/fo), we 
have l b 2^ — 1^/2. However, for K bi = iof exp(— nr 2 /l 2 ) we 

have l b 2 — lb i /y/2. There is no universal relation between the 
correlation scales of even these simply connected variables. Such a 
relation should be established in each specific case from the statis- 
tical properties of each physical component of the system. 

We calculate the correlation length l e of the synchrotron emis- 
sivity £ oc n C rB± in the synthetic radio maps from its autocorre- 
lation function K s (r), for various values of the cross-correlation 
coefficients C, Bo and no : 

h = a' 2 (£ K e (r) dr - < £ > 2 ) , (19) 

where L is the pathlength (assuming L 2> le) and o> is the standard 
deviation of the synchrotron emissivity. 

The resulting dependence of l e on the integral scale of the 
magnetic field l b = n/(2ko), obtained in Eq. J15b , can be approxi- 
mated as 

l e /L^G- x {Al b /L) - 6t , (20) 

where the numerical factor G depends on the model parameters 
and the cross-correlation coefficient C. The contours of G in the 
(S„, C)-plane are shown in Fig. [7} G — 9-10 are representative 
values for 5 b > 2-3. The resulting values of TV — L/(2l e ) are 
used below to compare the synchrotron fluctuations obtained from 
observations in Section|4]with the model of Section[5] 
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Figure 8. Relative fluctuations of the synchrotron intensity in synthetic ra- 
dio maps of Section |6~71 versus the cross-correlation coefficient between 
B 2 and n cr for various choices of the model parameters. Top panel: a se- 
lection of S b values for fixed S n = 0.5 (solid: S b = 10; dashed: 1; dotted: 
0.1). Grey horizontal lines con'espond to uniformly distributed cosmic rays, 
8 n = 0, for the same values of S b (with 8[ decreasing with 8 b ). Lower 
panel: various values of <5 n for S b = 10 (<5 n = 1, solid; 0.6, dashed; 0.4, 
dotted; 0.2, dash-dotted). 



6.3 Correlation between cosmic rays and magnetic fields 

The relative intensity of synchrotron fluctuations is sensitive to the 
number TV of correlations cells of synchrotron emissivity within the 
beam (or along the line of sight in case of a pencil beam). When 
comparing the theoretical model with observations, we adopt L — 
1 kpc for the pathlength in the Milky Way, l bi — 50 pc for the 
correlation length of magnetic field, 5t = 3 (the asymptotic limit 
S b 2> 1 is quite accurate in this case), and explore the range — 1 ^ 
C ^ 1 for the cross-correlation coefficient between cosmic-ray and 
magnetic fluctuations. For h/L — 0.05 and G ~ 9 (see Fig. O, 
we have l £ ~ 40 pc and N = L/(2l e ) ~ 10; we also discuss the 
effect of larger values of N. 

Figure [8] shows the dependence of 8i — 07 /To on the cross- 
correlation coefficient C for various values of the parameters S b 
and 8 n . The calculations are based on 100 realizations of B, so the 
statistical errors of the mean values shown are negligible. 

As can be seen from the upper panel of Fig. [8] the relative 
magnitude of synchrotron fluctuations, 8i\ // ~N > 0.7, obtained for 
Sb ^ 1 and S n = 0.5, is stronger than what is observed in the 
Milky Way, 5i</N = 0.3-0.6 assuming N = 10. If TV = 20, the 
conservative observational estimate Si — 0.1-0.2 translates into 
5!\fN = 0.4-0.9, implying C < -0.6 at the upper end of this 
range. Thus, 8 n < 0.5 seems to be justified, unless N is signifi- 
cantly larger than 10 or, otherwise, S b < 1 (which is highly implau- 
sible). Since the estimate Si = 0.1-0.2 has been obtained for high 
Galactic latitudes, the path length is unlikely to be much longer than 
1 kpc, and the correlation length of the synchrotron fluctuations can 
hardly be much shorter than about 50 pc. Thus, excluding the case 



of simultaneously large L and small l E , we conclude that the dis- 
tribution of cosmic ray electrons is unlikely to have any significant 
variations at scales of order 50-100 pc. 

The lower panel in Fig. [8] where a range of values of 8 n are 
used with 8 b = 10, suggests that any positive correlation between 
cosmic rays and magnetic fields is only compatible with the obser- 
vational estimate 8i\fN = 0.3-0.6 (for N = 10) if S n < 0.2. In 
fact, a upper limit SivN = 0.9 (for TV = 20) might be achieved 
for ricr = const. The values of 5ivN in this case are shown with 
grey horizontal lines in the upper panel: for example, 5/vTV = 1 
is compatible with 8 b = 10. However, the lower values of syn- 
chrotron fluctuations in the Milky Way in the range 5iVN = 0.3- 
0.9 for TV = 10-20 can be compatible with the presence of fluc- 
tuations in cosmic ray density mildly anticorrelated with those in 
magnetic field. It is difficult to be precise here, but S n < 0.2 and 
C < —0.5 seems to be an acceptable combination of parameters. 



7 PROPAGATION OF COSMIC RAYS AND 
EQUIPARTITION WITH MAGNETIC FIELDS 

To illustrate the relation between cosmic rays and magnetic fields, 
consider a simple model of cosmic ray propagation near a mag- 
netic flux tube. The number density of cosmic rays n cr (or their 
energy density e cr ) is assumed to obey diffusion equation with the 

( Parker 



sourc e Q and diffusivity D depending on magnetic field 



1 19691 ; iKuznetsov & Ptuskirj 1 19831 : [Schlickeise r & Lerchd 



19851) . 



Consider a steady state of a one-dimensional system with Q = 
const. Magnetic field is assumed to have a statistically uniform flue- 
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tuating component, b 2 — a 2 = const, whereas the mean field is 
confined to a Gaussian slab of half-thickness L symmetric with re- 
spect to x = 0: By = B exp[-~x 2 /(2d 2 )], B x = B z = 0. The 
cosmic ray diffusivity is assumed to depend on the relative strength 
of magnetic fluctuations, D — Doa 2 /B 2 . The resulting steady- 
state diffusion equation 

4- D 4-n C r + Q = 0, 

da; di 

can easily be solved with the boundary conditions 

-^r(°) =0 > "crllsl-foo = 0, 

to yield 

2D <t£ 

The total number (and energy) of cosmic rays remains finite despite 
the uniform distribution of their sources, Q = const, because D — > 
oo as \x\ — > oo in this illustrative model. 

This simple solution shows that, near a magnetic flux tube in a 
statistically homogeneous random magnetic field, cosmic rays con- 
centrate where the total magnetic field is stronger because their dif- 
fusivity is smaller there. In this example, the spatial distributions of 
cosmic rays and magnetic field are tightly correlated. 

Another type of arguments relating cosmic ray energy den- 
sity to parameters of the interstellar medium was suggested by 
IPadoan & Scaiol J2005h . If both the magnetic flux and the cos- 
mic ray flux are conserved, BS = const and n CI US = const 
(where B is the magnetic field strength and S is the area within 
a fluid contour, and U is the cosmic ray streaming velocity), one 
obtains n CI U/B = const, which yields n cr oc n 1//2 , given that 
U = Va oc Bn" 1 ' 2 , with n the gas number density and Va the 
Alfven speed. Thus, cosmic ray energy density is independent of 
magnetic field strength and scales with thermal gas density. This 
result relies on the fact that the streaming velocity of the cosmic 
rays is proportional to the Alfven speed. If, instead, U = V, with 
V the gas speed, we obtain n cr oc B from these arguments. No 
clear sc aling of e cr w i th ma gnetic field was observed in the simula- 
tions of lSnodin et al.l 12006) who use the gas velocity for U. 

Assumption of the detailed, point-wise (local) equipartition 
between cosmic rays and magnetic fields is dubious also because 
these two quantities have vastly different diffusivities, and there- 
fore cannot be similarly distributed in space. Magnetic filaments 
and sheets produced by the small-scale dyna mo in the diffuse warm 
gas can have scales as small as a few parsecs JShukurovl2 007). and 
the strength of this turbulent magnetic field can be about 5 /iG. The 
large-scale magnetic field varies over scales of order 1 kpc, consis- 
tent with the turbulent diffusivity of 10 26 cm 2 s _1 and time scale 
5 x 10 s yr. The diffusive length scale of cosmic rays, based on the 
diffusivity of D ~ 10 28 cm 2 s _1 and the confinement time in the 
disc, r ~ 10 7 yr, is about (2Dr) 1/ ^ 2 ~ 1 kpc, similar to that of the 
large-scale magnetic field. On these grounds, is not impossible that 
the energy densities of cosmic rays and the large-scale magnetic 
field vary at similar scales, but this would be very implausible for 
the total magnetic field. Then equipartition arguments may be bet- 
ter applicable to observations of external galaxies, where the linear 
resolution is not better than a few hundred parsecs, than to the case 
of the Milky Way. 



8 DISCUSSION 

The general picture emerging from our results is where cosmic 
rays and magnetic fields are slightly anticorrelated at the relatively 
small, sub-kiloparsec scales explored here (n cr = const is also a 
viable possibility). Such an anticorrelation can result from the sta- 
tistical pressure equilibrium in the ISM, where cosmic rays and 
magnetic fields make similar contributions to the total pressure. 
An additional effect leading to an anticorrelation is the increase in 
the synchrotron losses of relativistic electrons in stronger magnetic 
field. However, equipartition between cosmic rays and magnetic 
field c annot be excluded at larger scales of order and larger than 
1 kpc. iHoernes et all (l998) indirectly make a similar conclusion 
concerning loss of equipartition at small scales from their analysis 
of the radio-far-infrared correlation in M31. 

Since magnetic fields and cosmic rays have vastly different 
diffusivities, and therefore, must vary at very different scales, any 
strong correlation between them can hardly be expected at scales 
smaller than 1 kpc. Correlated (or rather anticorrelated) fluctuations 
can, however, arise from such secondary processes as the adjust- 
ment to pressure equilibrium, etc. 

Our arguments and conclusions are based on observations and 
modelling of synchrotron emission, a tracer of the electron compo- 
nent of cosmic rays. Thus, our conclusions strictly apply to the cos- 
mic ray electrons. However, the only significant difference between 
the behaviour of electrons and protons in this context is that the 
former experience by far stronger energy losses due to synchrotron 
emission and inverse Compton scattering off the relic microwave 
photons. The energy loss time scale ~ 4 x 10 8 yx{E/l GeV) -1 in 
magnetic field of 5 /xG in strength, for particles emitting at wave- 
lengths larger than 1 cm, is much longer than the confinement time 
10 7 yr, so the energy losses are negligible unless the local magnetic 
field is unusually strong. Therefore, we extend our conclusions de- 
rived from analysis of synchrotron fluctuations to cosmic rays as a 
whole. Moreover, energy losses can only make the distribution of 
the electrons more inhomogeneous than that of the protons, so that 
our conclusions are robust with respect to this caveat. 

Our results can significantly change the interpretation of high- 
resolution radio observations of the Milky Way and spiral galaxies. 
Present interpretations, aimed at estimating the strength and geom- 
etry of interstellar magnetic fields, rely heavily on the assumption 
of local equipartition between cosmic rays and magnetic fields, at 
a scale corresponding to the resolution of the observations. This 
assumption is acceptable if the resolution is not finer than the dif- 
fusion length of the cosmic rays, about 1 kpc. However, this as- 
sumption is questionable when applied to observations at higher 
resolution. We can suggest a different procedure to interpret such 
observations. The original radio maps should first be smoothed to 
the scale of the cosmic ray distribution, 1 kpc, where the equipar- 
tition assumption may apply, and the distribution of cosmic rays 
can be deduced from the smoothed data. After that, this distribu- 
tion of cosmic rays can be used to deduce magnetic field distribu- 
tion from the data at the original resolution. Since a larger part of 
the synchrotron fluctuations observed will be attributed to magnetic 
fields, it is clear that this procedure will result in a more inhomoge- 
neous magnetic field than that arising from the assumption of local 
equipartition so often used now. 
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APPENDIX A: THE CROSS-CORRELATION 
COEFFICIENT OF B 2 AND n CT 

Here we derive analytical expressions for the mean value and the 
standard deviation of synchrotron intensity for the model of par- 
tially correlated magnetic fields and cosmic rays introduced in Sec- 
tion [5] We assume that each Cartesian vectorial component of the 
random magnetic field b is a Gaussian random variable with zero 
mean value, b; = 0, and that the vector field is isotropic, 



and 



where i — x,y,z and bar denotes ensemble averaging. Then 

(2k)\ 



u2k _ K^LL 2k ,2k+l _ n 

The cross-correlation coefficient of n cr and B 2 is defined as 



(Al) 



(A2) 



C(B 2 ,n c 



n CI B 2 -nB 

(JnCT B 2 



Given Eq. J17t . it can be shown that 

n cr B 2 = a<j 2 B i + noB 2 and er 2 = a 2 a 2 B 2 + /3 2 a 2 ^ 2 - 

Then C(B 2 ,n CT ) is equal to a given value C, 

a 



a 2 a 2 + f3 2 a 2 



= c, 



D- 



provided 



C 



Vi -c 2 



n = and P — 

o B 2 a B 2 



For given n cr and B, the mean value and the standard devi- 
ation of the synchrotron intensity can be calculated using the fol- 
lowing relations, obtained using Eqs dAlt and dA2t : 



~D4~ d4 , 8 2 2,8 4 

B l = B 0± + -B 0± a b + -<j b , 



B 2 Bi = B% ± B^ + -Bl^alBl 



8 4 D 2 



n6 . 19 c 4 8 .80 J 4 56 6 



B 4 B ± = 2S ox Bo|| + l^B 0± a b B \\ + B 0± B ^ 
64 2 4 D 2 ,8 2 2 D 4 

+ -g-- t >0± Cr i J -D0|| + g- t >0± <T 6- t '0|| 

16 6 D 2 . 8 4 4 s 1 34. g 2 

109 4 4 , 104 a 6 56 8 
+ — i>oi^ + —^~B 0± a b + -g-ffe- 

From Eq. (|6}, we have I = L {e) L , where angular brackets 
denote spatial averaging, with subscript indicating the averaging 
region: L for the line of sight, S for the area and V for the volume 
averaging. This leads to 7q = {I) s = L {e) v = Nl £ e provided 
that N = L/(2l s ) > 1. Finally, we obtain, from I = Ln CT B\ 
ando-f = Lleie 2 -e 2 ): 

Nl 



4 



Io = -rT t 4 "^' 7 * + 6a? + 3 (4aS n ai + 3) B, 



Ox] 5 



(A3) 



af = ^|224Q 2 ^ag + 144q5„ (9aS n B% ± + l) at 

+ 18 [5lQ 2 (5 2 i3o X + 36a6 n B 2 ± - 4^ (C 2 - 1) + 2 
+ 12a 2 «5 2 (9Bt )± + 24a 2 b B? )± + 8a 6 4 ) B 2 {l a 2 b 
+ 108 [{aS n B 2 ± + l) 2 - 2«S 2 (C 2 - 1)] B 2 ± a 2 

-815 2 (C 2 -l)B f tx}. (A4) 
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